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1. Introduction 

Luminous matter almost certainly does not account for all matter in the 
Universe [1], and this matter deficit inspires both particle-physics and astro- 
physics speculation. Among the particle-physics solutions discussed, perhaps 
the most attractive idea is that stable weakly interacting massive particles 
(WIMPs) may make up the dark matter. Currently, the most promising candi- 
date WIMP is the neutralino [2] [3], a linear combination of the supersymmetric 
partners of the photon, Z°, and Higgs bosons. 

Recently, we have presented a complete calculation of the cross-section for 
neutralino annihilation to the two-gluon final state [4]. This one-loop process 
can be competitive with tree-level annihilation processes in large regions of the 
parameter space of the supersymmetric standard model. The main objective 
of Ref. [4] was to calculate the resulting dilution of the secondary-neutrino 
spectrum, which is of interest since such secondary neutrinos can provide a 
detectable signal for annihilation of neutralinos captured in the Sun or the Earth 
[5]. A large branching ratio to the two-gluon final state implies a decreased 
branching ratio for the annihilation to heavy-fermion final states, which produce 
the signal secondary-neutrino spectrum via weak decay. 

In Ref. [4], the possibility was discussed that the gluonic final state, though 
decreasing the detectability via secondary neutrinos, may open another avenue 
of detection — low-energy cosmic-ray antiprotons from WIMP annihilation in 
the galactic halo [6]. Hadronization of the final-state gluon jets will produce 
antiprotons in some fraction of the annihilation events. Although the "back- 
ground" spectrum of cosmic-ray antiprotons produced by interactions of pri- 
mary cosmic rays with interstellar matter is quite uncertain, the antiproton 
spectrum in standard cosmic-ray propagation models falls steeply at energies 
less than about a GeV [7]. If neutralinos populate the halo, the resulting distor- 
tion of the background spectrum due to antiprotons from WIMP annihilation 
may be significant at low energies. Previous work has considered the possibility 
that a detectable cosmic-ray antiproton signal could arise from the hadroniza- 
tion of b-quark and c-quark jets in the tree-level annihilation processes [8] [9], 
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and we must include this mode for antiproton production as well, though our 
discussion will concentrate on the new gluonic contribution. 

Due to numerous astrophysical uncertainties, it is impossible to make a pre- 
cise prediction for the antiproton spectrum for a given particle-physics model. 
However, we can show, for reasonable astrophysical assumptions, that a low- 
energy cosmic-ray antiproton signal from WIMP annihilation may indeed be 
observable. 

In the following we will give an estimate of the source spectrum for an- 
tiprotons produced by neutralino annihilation. Then we will present numerical 
results for the antiproton spectrum at the Earth, for some model parameters. 
We concentrate on a simple B-'mo WIMP scenario to illustrate the possibilities. 

2. Source Function for Cosmic-Ray Antiprotons 

In the annihilation of neutralinos to two gluons, XX ~~ * 99 1 the back-to-back 
gluonic jets, each with energy m^, will shower to produce color-singlet hadronic 
states, notably pions, kaons, (anti)protons, and (anti)neutrons. This process is 
described by fragmentation functions, D 9 ^ (z), one for each species of hadron, 
h, which give the distribution in energy of the final-state hadrons as a fraction 
of the initial jet energy, z = -E^ad/ -Ej e t- There is also a small logarithmic 
dependence on the jet energy, calculable in perturbative QCD [10]. Clearly the 
spectrum of injected antiprotons is precisely (twice) the fragmentation spectrum 
for the hadronization of a gluon jet, g — > pp + X. Therefore we must assemble 
what is known about gluonic jet fragmentation to (anti)protons. Note that we 
are interested in the properties of gluon jets with energies below about 150 GeV, 
since the branching ratio to the gluonic final state can be appreciable only when 
the neutralino mass is smaller than the top-quark mass and, in some cases, the 
W- and Z-boson masses. 

The fragmentation process is an inherently non-perturbative QCD process, 
beyond our ability to calculate at the present time. Therefore, we rely on the 
measured properties of QCD jets and on fragmentation models when necessary. 
Unfortunately, very few properties of gluon jets are measured well. Though 
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gluon jets are produced copiously in pp collisions, the required exclusive analysis 
of the final state in such collisions is extremely difficult. Many properties of 
quark jets have been measured in e + e~ collisions [11] [12] [13], but the analysis 
of three-jet events, attributed to the final state qqg, has not produced the 
required detailed gluon fragmentation information. However, a combination 
of theoretical arguments and experimental information can be used to relate 
the properties of gluon and quark jets. Furthermore, detailed Monte-Carlo 
fragmentation models provide some insight into gluonic jets; the reliability of 
such models is judged based on their ability to predict many aspects of jet 
physics in e + e~ collisions [12]. 

As an empirical fact, the fragmentation functions for jets produced in e + e~ 
annihilation can be parametrized in the form 

1 d<j(e + e~ — > ppX) v—, 

v — yy_L = j3c Vojexp-fciz, (2.1) 

i 

where C is a constant dependent on the jet energy which fixes the average 
multiplicity, and j3 is the velocity of the antiproton [12] [13]. Interpreted as 
the differential probability that a quark-jet pair will hadronize to produce a pp 
pair in the final state, such an expression is precisely (twice) the fragmenta- 
tion function for q — > p + X . Unfortunately, we need such a parameterization 
for gluon jets, and data is not available. However, there is some theoretical 
work which can help to relate the fragmentation of quark and gluon jets. It 
can be shown that gluon jets should have higher average multiplicities, in the 
asymptotic energy regime [14]. Also, successful Monte Carlo models confirm 
the expectation of higher average multiplicities for gluon jets [15]. There is 
also some expectation that gluon jets give rise to softer hadron spectra [10] [16] . 
However, this effect is not large, and we will choose the fragmentation functions 
for gluon jets to have the same form as for quark jets, corrected only for a larger 
average multiplicity by choice of a different normalization constant. The effect 
of this prescription will be a slight underestimation of the low-energy antipro- 
ton spectrum, and this is conservative for our purposes. Furthermore, some 
fragmentation models can successfully reproduce the observations assuming no 
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difference between gluon- and quark-jet fragmentation [17], and thus we do not 
consider this prescription to be the leading source of error in our estimate. 

Our fit to the data for e + e~ — > ppX, as reproduced in Ref. [12] and Ref. 
[15], gives a quark-jet fragmentation function of the form 

D p (z) = C(3 [exp(-25z) + 0.1 exp(-5.3z)] . (2.2) 

Note that this data includes contamination from the untagged b — > pX frag- 
mentation processes which are expected to have a somewhat different functional 
form; however, this effect is small. 

In order to fix the constant C, and thus normalize the gluon fragmenta- 
tion function, we make use of the fact that the average multiplicity of charged 
particles in jets follows the approximate formula [12] [18] 



(n ch ) = a + &exp 



.1 S ^ 

c[ln Ql 



(2.3) 



We find that a ~ 2.8, b ~ 0.05, c ~ 2.2, and Q ~ 1 GeV provides an 
adequate description of the multiplicity in the Webber model, as illustrated in 
Figure 10 of Ref. [15], in the jet energy regime in which we are interested. This 
gives the "per event" multiplicity rather than the "per jet" multiplicity, and 
thus the factor of two which counts the number of jets is implicit. In addition, 
we have included the jet multiplicity data for jet events on the Z° peak [19] [11]. 

Data on exclusive quark-jet multiplicities [20], reproduced in Figure 9 of 
Ref. [15], together with the Webber model, indicate that the pp multiplicity is 
an essentially constant fraction of the total charged-particle multiplicity over 
the jet energy range of interest to us, (n p p) ~ 0.06(n c h), with an accuracy of 
about twenty-five percent. Note also that the production of antineutrons which 
decay in flight gives a factor of two increase, which we choose to include in the 
normalization constant rather than display explicitly. 

Finally, we must include the contributions to the antiproton spectrum from 
fragmentation of the bb and cc final states of neutralino annihilation, which have 
been discussed previously in the literature [8] [9]. We will use the b, c — > pX 
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fragmentation functions, Dp q (z), of Ref. [9]. The factor of two due to anti- 
neutron production is included in the normalization of these functions. 

The gluon fragmentation function has been normalized by J Dp 9 (z) dz = 
(n p p) g , and the normalization of the heavy quark fragmentation functions is 
included in the parameterization of Ref. [9]. As discussed in Ref. [9], the nor- 
malization of the fragmentation functions must be evolved for jet energies away 
from 10 GeV, where the fit parameters are quoted; this evolution is numerically 
significant. 

Defining the antiproton source function, qp, as the differential number of 
produced antiprotons per unit volume, per time interval, per energy interval, 
we write 

qp(E) = q$^"(E) + qf^ q ~(E), (2.4) 

where qq indicates the heavy quark final states, bb and cc. At the small velocities 
typical of WIMPs in the halo, annihilation of neutralinos to light quarks is 
negligible. It is a simple matter to write the source functions for antiprotons 
produced in neutralino annihilation in terms of the fragmentation functions and 
the annihilation cross sections 

q^"(E p ) = (av) xx ^ gg {n\) ±- D 99 (z), (2.5) 

q^"(E p ) = (av) xx ^ qq - {n\) ±- D*\z), (2.6) 

m x 

where z = Ep/m x is the fraction of the jet energy carried by the antiproton, and 
is the volume average of the square of the number density of neutralinos 
in the galactic halo. 

Locally, ?i^ >cal = 0.4 (m^/GeV) _1 po.4 cm -3 , where po.4 is the local halo 
density in units of 0.4 GeV cm -3 . However, if there is a bulge population of 
WIMPs and the halo core radius r COTe is small compared with R ~ 8 — 10 kpc, 
our distance from the center of the galaxy [21], then there can be a significant 
enhancement in the volume average of the square of the number density [22], 

(4) - T^^ OCaI ) 2 ' forfl » r - 

J-^' core 



Now, for illustration, we consider the simplified B-'vuo model discussed in 
Ref. [4]. If the neutralino is primarily B-'vuo and is less massive than the top 
quark, then annihilation into gluons can be significant. In addition, there are 
some theoretical arguments that suggest that the neutralino is likely to be a 
B-'vuo in this mass range [23]. The neutralino relic abundance Q x (in units of 
the critical density) is determined by the cross section for annihilation of two 
neutralinos into all lighter particles. The thermally averaged cross section for 
annihilation of two neutralinos, at a temperature T in the early Universe, is 
generally written 

<™>ann-°+ 6 ( T A"*)- (2-7) 

Here, a is the s-wave contribution to av and is the only part which survives in 
the limit v — > 0; b is the p-wave contribution. Neutralinos "freeze out" when 
T/m x ~ 1/20. If a <C b(T/m x ) when the neutralinos freeze out, then the relic 
abundance is roughly (see, e.g., Ref. [24]) 

Q x h 2 ~ ( h - j\ , (2.8) 

x \8.8 x 10" 9 GeV -2 / V ' 

where h is the Hubble parameter in units of 100 km sec -1 Mpc -1 . 

If the neutralino is a pure B-'mo with mass large compared to the masses 
of the light fermions, but small compared to the masses of the top quark and 
the sfermions, mj <C ra x <C m-t,mq, then the cross-section coefficients are 

, ( 100 GevV 2 , N 

o~4.9x 10" 7 GeV" 2 , (2.9) 

V m f J 

aI1C ^ 2 2 

fi f m x \ { 100 GeV\ 2 
6 ~ 1.2 x 10 -6 — GeV -2 . (2.10) 

\ m fJ V m f J 

Because of the helicity suppression of majorana-fermion annihilation into light 
fermions [25], a is suppressed by the squares of the quark and lepton masses, 
so a <C b. Then, the cosmological abundance of a pure B-'mo in the given mass 
range is roughly 



This implies a relation between the B-'mo mass and (common) sfermion 
mass, m x ~ 0.2 (fi^/i 2 /0.25) -1 / 2 (ra^/100 GeV). Neutralinos can account 
for the dark matter in galactic halos if the neutralino relic abundance falls in 
the range 0.03 < fi^/i 2 < 1; if neutralinos are to contribute closure density, 
then fi^/i 2 ~ 0.25. 

At this point we have specified all the parameters necessary to calculate 
the source spectrum of antiprotons. Recall the approximate cross-sections for 
the various channels in the pure B-'mo limit, as discussed in Ref. [4]. Because 
the two-gluon cross-section depends on the same combination of dimensionful 
parameters as the relic density, the B-'mo relic abundance actually fixes the 
two-gluon cross- section, a gg v ~ 4 x 10" 12 (Q^h 2 /0.25)- 1 GeV" 2 . Also, the 
cross-sections for annihilation to heavy quarks in the non-relativistic limit can 
be written as 

°» 

a 99 

&CC 

a 99 

Notice that the fraction of annihilations to quarks is relatively small compared 
to that to all fermions, since the latter is generally dominated by annihilation 
to tt. In addition, in (2.12), we have neglected the effect of running quark 
masses which would further decrease the contribution of the heavy-quark final 
states. Thus the gluon final state can be quite important, in comparison to the 
quark-antiquark final states. 

An interesting feature of the gluon contribution to the antiproton pro- 
duction is that the m x dependence of the cross-section cancels that which is 
implicit in the denominator due to the factors of the neutralino number density. 
In other words, if we assume, as we have done, that the neutralinos make up the 
dark matter in the galactic halo, with (local) mass density p ~ 0.4 GeV/cm 3 , 
then their number density in the halo is inversely proportional to their mass. 
This factor of m x 2 in the denominator is canceled by a similar factor in the 
two-gluon cross-section. On the other hand, there is no strong m x dependence 
in the cross section for annihilation into quarks, so the gluon contribution does 



30 GeV \ 2 

m x ) 
25 GeV\ 2 

m x ) 



(2.12) 
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not decrease as quickly with increasing neutralino mass. Therefore the gluon 
contribution is relatively more important for larger neutralino masses. Also, 
there are several other arguments for a relatively important gluon contribution. 
First, the average charged multiplicity for gluon jets is measurably higher than 
that for quark jets. Second is a more subtle effect; as an empirical fact, the frag- 
mentation production of a heavy-flavor meson in a heavy-flavor jet is relatively 
hard [12], compared to the similar process in a light quark jet. On average, 
this means that less energy is available for baryon production in heavy quark 
jets. On the contrary, a gluon jet fragments in much the same way as a light 
quark jet. Clearly the combination of these three effects can be numerically 
significant for the resulting cosmic-ray antiproton spectrum. 

3. Results for the Observed Spectrum 

Given the above source function q p (E) for antiprotons injected in the 
galaxy, the spectrum at the boundary of the heliosphere, <p^(E), is determined 
by a model for cosmic-ray propagation through the interstellar medium. The 
simplest model for interstellar propagation, which we choose, is the leaky-box 
model [7]. (An alternative which would increase the flux of background an- 
tiprotons, the closed-galaxy model [26], seems to be in conflict with recent 
measurements of the low-energy cosmic-ray 3 He/ 4 He ratio [27].) The distortion 
of the particle spectrum due to propagation in the leaky-box model is described 
by a (possibly energy dependent) escape time, T esc (E p ), and we write 

<l>f{E p ) = ±-cT esc {E p )q p {E p ). (3.1) 

The diffusion length contains the effects of galactic escape as well as depletion of 
the spectrum by interactions, which can be important for less energetic particles. 
We use r esc = 10 8 yrs [28], although we note that estimates for this value vary 
by roughly an order of magnitude. 
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To obtain the observed spectrum, <f>p bs , the interstellar spectrum must be 
further corrected for solar modulation. Following Ref. [9], we use the solar mod- 
ulation model of Perko [29]. The model of Ref. [29] transforms the interstellar 
spectrum to the observed spectrum by 



4>f s {E) = 4>f{E'), (3.2) 

where [9] [29] 



^^'-{pJn^^ + E^AE, for p < p c . 



(3.3) 



Here p c = 1.015 GeV, p = (E 2 — m 2 ) 1 ! 2 , where m v is the proton mass, and AE 
is an energy shift which depends on the phase in the solar cycle. 

The results for the modulated cosmic-ray antiproton spectra at solar min- 
imum, expressed as a fraction of the observed cosmic-ray proton spectrum, are 
displayed in Fig. 1. We took the interstellar cosmic-ray proton spectrum to 
be <f> v = 1.93 x 10 4 P(E/ GeV) -2 - 7 m" 2 s" 1 sr" 1 GeV" 1 [31], and modulated 
it in the same way as the antiprotons. The upper (lower) solid curve is the 
spectrum expected if a 30-GeV (60-GeV) B-'mo populates the halo. Here we 
have made the simplest assumptions; we have taken the relic density of WIMPs 
to be unity, and we have assumed (n^) = (n}~ cal ) 2 . The upper (lower) dashed 
curve is that which would have been obtained by ignoring the contribution 
from the two-gluon final state for a 30-GeV (60-GeV) B-'mo. The dotted line 
is the "background" flux of cosmic-ray antiprotons expected from spallation of 
primary cosmic-rays in the leaky-box model [26] suitably modulated at solar 
minimum. Also shown is the current observational limit of p/p < 2 x 10 -5 
in the kinetic-energy range 100-1500 MeV [30]. As seen, the two-gluon final 
state can be the dominant source of antiprotons if the neutralino is primarily 
B-'mo. In addition, Fig. 1 illustrates that WIMP annihilation could plausibly 
provide a low-energy cosmic-ray antiproton signal larger than that expected 
from background, yet small enough to have evaded detection thus far. 
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Fig. 1 Observed antiproton/proton ratio as a function of kinetic energy. The 
dotted curve is the expected solar-modulated flux of "background" antiprotons 
from standard propagation. The upper (lower) solid curve is the spectrum 
predicted if a 30-GeV (60-GeV) B-'mo populates the galactic halo. The upper 
(lower) dashed curve is the spectrum that would be predicted for a 30-GeV 
(60-GeV) B-'mo if we had ignored the contribution of the two-gluon final state. 
Also shown are the current observational upper limits [30]. 

4. Discussion 

The estimate presented here contains some uncertainties from the physics of 
hadronization as well as significant uncertainties from the physics of low-energy 
cosmic-ray antiproton propagation. By considering the source spectrum in some 
detail we hope that we have illustrated the nature of our results divorced from 
the latter uncertainties. It would seem difficult to use non-observation of such 
low-energy antiprotons to constrain dark-matter candidates. Unfortunately, 
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detailed improvements to the calculations are not easily implemented. For in- 
stance, many alternate models of cosmic-ray propagation can result in enhance- 
ments of the low-energy antiproton flux that might resemble the proposed signal 
from WIMP annihilation [7] [26]. Effects attributable to these models must be 
separated from any possible new-physics contribution such as neutralino anni- 
hilation. Hopefully, improved cosmic-ray data will help discriminate between 
the various models for secondary antiprotons. 

We should also point out that although the WIMP-induced cosmic-ray 
antiproton spectrum might indeed be observable, such a signal is by no means 
generic to models where WIMPS account for the dark matter in our halo. 
There is a relatively limited range of masses and cross sections for which the 
resulting relic density is large enough and for which the antiproton production 
cross section is high enough to produce an observable effect. Contrariwise, 
there are a number of reasonable alterations to the models considered here 
that could lead to a larger antiproton flux. Once again, we emphasize that 
the antiproton signal is likely to be enhanced for models where the energetic- 
neutrino signal from WIMP annihilation in the Sun or Earth is likely to be 
depleted. Therefore, the search for low-energy cosmic-ray antiprotons provides 
an important complement to energetic-neutrino searches in the effort to infer 
the existence of particle dark matter. 
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